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Abstract. Many more nuclear transitions have to be known in the determination of stellar reactivities for trans-iron nucleosynthesis 
than for reactions of light nuclei. This requires different theoretical and experimental approaches. Some of the issues specific for 
trans-iron nucleosynthesis are discussed. 


INTRODUCTION 

Nuclei beyond Fe are thought to have been created in a number of processes and in several different sites. The bulk of 
these intermediate mass and heavy nuclides were made in neutron capture processes, avoiding the problem of having 
to overcome high Coulomb barriers which otherwise suppress nuclear reactions in stellar plasmas. Historically, two 
main processes have been distinguished, the s-process (slow neutron capture) and the r-process (rapid neutron capture), 
each contributing about half of the abundances found beyond Fe. In the current understanding, however, these are split 
into subprocesses occuring in different sites [1]. For the s-process, there is a main component made in AGB stars and 
a weak component from massive stars [2]. The main component itself is made in two distinct environments inside 
the AGB star, during He-shell flashes and in interburst, intershell H-burning. The origin of r-process nuclei is less 
well understood but there are indications that also contributions from at least two different sites may be required to 
reproduce the r-process patterns found in the Sun and in stars of different ages [3, 4]. While the long-time favored 
core-collapse supernovae (ccSN) were found to have problems producing r-process nuclei, outflows from neutron-star 
mergers and magnetically driven jets from supernova explosions have attracted recent attention because models show 
promise for obtaining the required conditions for an r-process. A small number (32 - 35, depending on the model 
used) of proton-rich nuclei cannot be made in either s- and r-process and require other production mechanisms [5]. 
For most of these it has been shown that they can be made by photodisintegration of pre-existing nuclei in the outer 
layers of a massive star before and during its ccSN explosion. The lightest of these so-called p-nuclei, however, cannot 
be produced in sufficient quantities in this y-process and call for another explanation. Various suggestions have been 
made, such as a y-process in type la supernovae (thermonuclear disruption of a White Dwarf in a binary system), 
a vp-process in the deep layers of a ccSN, or a rp-process (rapid proton capture) on the surface of mass-accreting 
neutron stars. Finally, ^^^La (and perhaps partially also '*®Ta) has been shown to be produced by v-induced decays, 
again in the outer layers of a massive star during its explosion. 

Common to all these processes is that the required conditions (temperature range, matter density, neutron-to- 
proton ratio, duration) are constrained by observed abundances and by nuclear physics (reactivities, which are in 
turn determined by masses. Coulomb barriers, spectroscopy, low-energy y- and particle strength functions). Thus, it 
makes sense to classify nucleosynthesis by processes, for which the appropriate conditions may be found in more 
than one astrophysical site. This equally applies, of course, to stellar burning stages involving light nuclei. There 
are several important differences for heavier nuclei, however, which also affect which nuclear properties are relevant 
and what experimental and theoretical approaches can be used to determine astrophysical reaction rates. Special 
for nucleosynthesis beyond Fe are (i) generally higher plasma temperatures, (ii) higher matter density during the 
astrophysical burning process, (iii) higher intrinsic nuclear level density of the involved intermediate and heavy nuclei, 
and (iv) higher Coulomb barriers. 



Items (i) and (iv) result in higher effective interaction energies than in light element nucleosynthesis (but still 
low by nuclear physics standards), ranging from 8-30 keV for the main s-process and up to 90 keV for the weak 
component. Neutron captures in the r-process occur at around 80-120 keV. In charged particle reactions the effective 
energies are shifted to a few MeV, depending on the specific process and the charge of the nucleus and projectile (up 
to ss 9 MeV for a-reactions in the production of p-nuclei). Points (i) and (ii) also imply that (except for the s-process) 
short-lived nuclides are involved, which can be studied only in a limited manner in the laboratory, if at all. The higher 
inherent nuclear level density of heavier nuclei combined with the higher temperature leads - items (i) and (iii) - to 
considerable contributions of transitions from excited states of the target nuclei (see next section). 

Put in simple terms, these circumstances conspire to simplify a theoretical treatment (with exception of reactions 
on magic nuclei and close to the driplines) and complicate an experimental constraint of stellar rates. The many 
transitions (from target states to final states, mostly via compound states) involved in explosive nucleosynthesis lend 
themselves to the use of averaged quantities in their prediction. On the other hand, the large number of transitions 
restricts the applicability of direct and indirect experimental approaches studying a few transitions, as usually applied 
in the study of light nuclei, even when dealing with stable nuclides. Finally, only reactions on targets in the ground 
state (g.s.) are accessible in the laboratory whereas heavy element nucleosynthesis not only involves highly unstable 
targets but also such in excited states. 

In the following, some details of the above issues are briefly explained. 


Stellar cross sections and importance of thermally excited nuclei 


The stellar reaction rate (in reactions per time) for a reaction a + A —> b + B 'm?L plasma of temperature T can be 
expressed as [6, 7, 8] 
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with the number densities ria, tia, the reduced mass and Boltzmann constant k. The above equation also defines 
the reactivity R* - {o-*v)uA^bB which is often used synonymously with the term “rate”. Unfamiliarly in the nuclear 
physics context, the stellar cross section cr*(£, T) does not only depend on interaction energy but also on the plasma 
temperature which determines the population of excited states in the target. Explicitly showing the transitions between 
states i in nucleus A and final states j in B, it can be written as [6] 
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The partial cross sections cr'^^ are evaluated at the energy E - E,, noting the convention by [7] to zero cross sections 
at non-positive energy. Target levels are characterized by their spin 7, and excitation energy with the g.s. having 
i - 0 and £o = 0 MeV. The nuclear partition function of the target is 

G{T) = ^ (27,- + 1) ^ P,- . (3) 


It is often overlooked that the relative importance of transitions from excited states is given by the weights 
Wi - \ - EJE, which show a linear dependence on excitation energy, despite of the exponential decline of the thermal 
Boltzmann population factors P,. Relying solely on the P, as indication of the importance of thermally excited states 
will strongly underestimate their contribution. 

Another instructive way to judge the importance of excited target state transitions is to look at the g.s. contribution 
to the stellar rate which is defined as [10] 
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with cro(£) = Yjj being the g.s. reaction cross section as usually determined experimentally. It always holds that 
0 < Vo < 1. It is very important to note that this is different from the simple ratio Rq/R* of g.s. and stellar reactivities. 




FIGURE 1. Ground state contribution X to stellar (n, 7 ) rates at 0.384 GK (left) and at 1 GK (right) [9]. 


respectively, which has been used mistakenly in the past to quantify excited state contributions. As an example. Fig. 1 
shows Xq for neutron captures at 0.384 GK (30 keV effective energy, typical for the s-process) and 1 GK (typical for 
the r-process). It can be seen easily how the importance of excited state transitions increases for heavier nuclei. Even 
at s-process temperatures g.s. contributions can already be small, especially in the rare-earth region. The implications 
of this concerning experimentally unconstrained rate uncertainties have been discussed in [10]. The g.s. contributions 
also decrease strongly with increasing temperature. The only exceptions are reactions close to magic numbers which 
retain larger g.s. contributions also at higher T owing to the larger level spacing. Regarding charged particle reactions, 
Xq are even smaller in most cases, due to the shift of effective interaction energy to higher energies in the integration 
for the rate (Eq. 1) [11]. All Xq are given in [9]. 

As a rule of thumb when comparing g.s. contributions between forward reaction a + A —» b + B and reverse 
reaction b + B ^ a + A, the reaction direction with positive reaction Q-value almost always shows larger Xq. Among 
the comparatively few exceptions (due to the Coulomb suppression effect of excited state transitions [12]) are capture 
reactions, which always have much larger Xq than their photodisintegration counterparts. In fact, Xq < 5 x 10“^ for 
astrophysically relevant photodisintegration rates [6, 9] (see also Section III.C.l in [1] for a detailed discussion of the 
problematic application of photodisintegration experiments to astrophysical topics). 

The Xq can also be used to calculate a new uncertainty when combining experimental and theoretical information 
in an improved stellar rate and thus demonstrate how much a stellar rate can be constrained by a measurement at best 
[1]. Refs. [9, 10] also provide the recommended procedure for including measured ctq in an improved stellar rate. In 
a nutshell, instead of correcting the experimental reactivity by some factor accounting for the thermal population of 
excited states, the theoretical stellar reactivity B* has to be corrected by a factor /* accounting for the experimen¬ 
tal information on the g.s. reactivity The new stellar reactivity combining theory and 

experiment, is obtained by using the factor 
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Obviously, also the remaining uncertainty of the stellar reactivity is affected by this procedure. The new uncertainty 
factor of is constructed from the original (theory) uncertainty factor m* > 1 and the experimental uncertainty 
factor t/exp > 1 of applying 


«:ew(7’) = t/exp + (m*(T) - t/exp)(l - Xo(r)) 
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Reaction mechanisms and important nuclear properties 

Apart from increasing the contribution of excited state transitions in the stellar rate, the fact that heavier nuclei have a 
higher nuclear level density (ED) and thus smaller level spacing than light species has another important consequence. 


























The LD at the formation energy ftform of the compound nucleus a+A-C-b + B determines the dominant reaction 
mechanism. In the absence of levels to be populated close to the given energy, direct reactions to low-lying final 
states dominate the reaction cross section [6, 14, 15]. When there is a small number of well separated excited states 
close to £form in C, this gives rise to resonances in the reaction cross sections, which can be described by the Breit- 
Wigner formula or by the R-matrix method [6, 16]. The considerable challenge therein is to determine the properties 
of the resonances contributing to the reaction rate integral. It is experimentally challenging to perform the required 
measurements for unstable nuclei and/or at low energy and theoretical ab initio methods cannot be applied to heavy 
nuclei, yet. 

Finally, an extreme case of resonant reactions appears when the LD is high, leading to a large number of un¬ 
resolved resonances which can be described by averaged resonance parameters. This is called the Hauser-Feshbach 
approach [17, 18]. The vast majority of reactions in nucleosynthesis of intermediate and heavy nuclei can be de¬ 
scribed in this model [6, 19] Despite of the increased number of nuclear transitions to be included in the model, this 
facilitates the predictions somewhat, as averaged quantities can be used. Calculations of reaction rates from smooth 
Hauser-Feshbach cross sections are also somewhat more “forgiving” to fluctuations around the “true” cross section 
value because of the integration over the projectile energy distribution. Different quantities, though, have to be known 
as input to the calculations than for direct or resonant reactions, such as optical potentials (related to the effective 
interaction in a many-nucleon system) and level densities (but not, e.g., spectroscopic factors of isolated levels) [6]. 
Nevertheless, Breit-Wigner and Hauser-Feshbach cross sections are closely related and the expressions for the cross 
sections look very similar [6]. Also, the reactivities depend on similar expressions in both cases, a fraction 
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The quantities W'^, W^, and are the total width of a compound level in C with spin Jc at excitation energy Ec 
and including all transitions to levels i in nucleus A with £, < Ec - 5 pro in the entrance channel, the total width of 
this compound state including all possible transitions to levels j in B with Ej < Ec - Sin the exit channel, and the 
combined width = Tja^a.b,... ^Ec Jc emission processes from the compound state to all 

channels A, respectively [8, 6]. The separation energies of the projectile and the ejectile in the compound nucleus C 
are denoted by 5pro and S^j, respectively. Assuming narrow resonances, the rate of a single Breit-Wigner resonance at 
excitation energy Ec (this would be at an energy of E = Ec - 5 pro in the entrance channel) is given by [6, 9] 
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In the Hauser-Feshbach model it is assumed that there is a large number of resonances at any compound formation 
energy Ec and for any Jc and parity ttc- Therefore, the rate is 

^HF ^ f J^r(E + Sp,o,Jc,trc)Ee-^^(‘^^^dE . (9) 


In this case, the widths used in 'F are assumed to be averaged quantities for a large number of similar resonances. 

The similarity of F as defined in Eq. (7) in the Breit-Wigner and Hauser-Feshbach approaches leads to the 
expectation of similar dependences on input quantities, i.e., on (averaged) widths. This is actually the case and allows 
to better disentangle the seemingly complicated sensitivity of Hauser-Feshbach predictions on input quantities such as 
masses, optical potentials, y-strength functions, and LDs. Often, these dependences are explored by arbitrarily using 
a selection of predefined treatments of these quantities and observing the change in the resulting cross sections or 
reactivities. This is then taken as an indication of the uncertainty in the underlying model. From a theory point of 
view, this is an inadequate or even incorrect procedure as an arbitrary selection of input models can never catch the 
actual uncertainty (see [9] for error estimates in predictions). An uneducated use of models may be misleading also 
because some of the models may have known limitations, discouraging the application to the chosen reaction. Rather, 
a systematic variation of the fundamental quantities, i.e., the reaction widths, is called for. Similarly, when trying to 
model a measured cross section, the starting point should be systematic variations of widths. After having identified 
the actually relevant widths and their required changes, it can be discussed how to obtain these based on theory. 

For systematic sensitivity studies it is convenient to define the sensitivity [9] 
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FIGURE 2. The maximally contributing El y energies when capturing 60 keV neutrons on Sn and Pb isotopes with even mass 
numbers are compared to the neutron separation energies 5n in the compound nuclei [13]. Note that the horizontal axis gives the 
neutron number N of the compound (final) nucleus. 


It is a measure of a change by a factor of vn = Qnew/^^oid in the output Q. as the result of a change in the input quantity 
q by the factor Vq - ^new/^oid, with s = 0 when no change occurs and s = 1 when the final result changes by the same 
factor as used in the variation of q, i.e., s = 1 implies vn = v^. Negative s are usually taken to indicate that q and Q vary 
in opposite directions. A large scale study of sensitivities for reactions with light projectiles from proton- to neutron- 
dripline has been published in [9], where g.s. cross sections and stellar reactivities have been chosen as Q and particle- 
and y-widths W, as well as LDs, were varied as q. While it may not be surprising that neutron captures are mainly 
sensitive to variations in the y-widths, due to the y-width being much smaller than the neutron width, this does not 
apply to reactions with charged particles. In stellar reactivities, charged particle widths are always smaller than both 
neutron- and y-widths, due to the small interaction energies and high Coulomb barriers. This poses a big challenge 
for experiments as it is very difficult to impossible to determine charged particle widths at astrophysically relevant 
energies. Another finding was that only those reactivities are sensitive to changes in the LD which are also sensitive 
to the y-width. This is because transitions to low-lying, discrete excited states dominate the particle widths (due to 
their larger relative energies) whereas y-widths are dominated by transitions with smaller relative energy to levels with 
higher excitation energies [6]. The latter can be understood by realizing that the competition between increasing LD 
with increasing excitation energy in C and decreasing y-strength (because of decreasing relative y-energy in C) leads 
to a peak in the contribution of single y-transitions at a specific energy Interestingly, ss 2 - 4 MeV for 

astrophysical neutron- as well as charged particle capture [13]. This rule is only broken for reactions with low isform 
in nuclei with low LD, such as at magic numbers or close to driplines, where y-transitions to the g.s. dominate. Figure 
2 shows an example for neutron captures on even Sn and Pb isotopes. More details are found in [13]. The knowledge 
of the contributing y-energies is important for judging the relevance of modifications to the photon-strength function 
(such as pygmy resonances) and the LD. Changes will affect reactivities only when they appear at energies around 
E'^^, otherwise they will be inconsequential for astrophysical applications. 

Exceptions to the applicability of the Hauser-Feshbach model are reactions with very low or negative projectile 
separation energy in C because then Eform is shifted to very low excitation energies, at which the LD may be too low to 
apply the model even for heavier nuclei. Typically this occurs close to the driplines. Moreover, due to the lower LD at 
magic nucleon numbers, the Hauser-Feshbach model also may not be applicable to magic target nuclei and low plasma 
temperatures [19]. This applies mainly for neutron captures, however, because the projectile energy range relevant for 
the calculation of the reaction rate is shifted to higher compound excitation energies when there are charged particles 











in entrance or exit channel [11]. 

The importance of direct reactions - and specifically of direct neutron capture - in astrophysical environments 
has been discussed in literature (see, e.g., [6, 14, 15] and references therein). Their sensitivities, however, are compar¬ 
atively simple to estimate because only single transitions are involved instead of summed widths. Thus, the rate of a 
direct reaction will linearly depend on the spectroscopic factors used. They will also be more strongly sensitive to a 
change in the properties of the single final state (excitation energy, spin, parity, and whether it is bound or unbound) 
than compound reactions [15]. Direct reactions at high energies - (d,p) and (d,n) reactions at several tens to hun¬ 
dreds of MeV, for instance - have also been suggested to be used in nuclear spectroscopy studies of unstable nuclei 
at radioactive ion beam facilities. Cross sections of reactions at such high energies do not appear in the calculation 
of astrophysical reaction rates but the extracted information on low-lying excited states and spectroscopic factors is 
useful in the calculation of the widths W implicitly appearing in Equations (8) and (9). 


Conclusion 

Many more transitions have to be known to determine stellar reactivities for trans-iron nucleosynthesis than in light 
element nucleosynthesis. This requires different theoretical and experimental approaches. Due to pronounced thermal 
plasma effects in the stellar rates, most measurements can only support predictions by testing models for single 
transitions but not constrain a stellar rate independently. It has to be made sure that astrophysically relevant transitions 
are studied, though. Published systematic sensitivity studies and g.s. contributions to the stellar rate help to guide 
experiments. On the other hand, predictions are simplified by being able to average over many transitions and apply 
the Hauser-Feshbach model, which has been successful in describing a large number of reaction cross sections. The 
challenges for theory are twofold. First, nuclear structure models have to be improved to be able to reliably predict the 
nuclear properties (such as nuclear spectroscopy, FD, or y-strength functions) required for the astrophysical reaction 
cross section prediction. Another major challenge to theory is to accurately describe the competition between direct, 
resonant, and Hauser-Feshbach reaction mechanisms for nuclei close to magic numbers and close to the driplines. 
Some first attempts have been made to combine direct and Hauser-Feshbach cross sections for neutron-rich nuclei 
[6, 20] but currently a predictive treatment of individual resonances (which nevertheless may be very important also 
for magic nuclei and close to driplines) is beyond the reach of theory. 
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